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Abstract 

The equations governing the interaction between incoherent light and electron-acoustic waves 
are presented. The modulational instability properties of the system are studied, and the effect 
of partially coherent light is discussed. It is shown that partial coherence of the light suppresses 
the modulational instability. However, short wavelength perturbations are less affected, and will 
therefore dominate in, e.g. pulse filamentation. The results may be of importance to space plasmas 
and laser-plasma systems. 
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The investigation of the interaction between light and ion-acoustic waves was pioneered by 
Karpman in the early seventies □ D . Since then there have been numerous studies devoted 
to the corresponding governing equations, which are known as the Karpman equations^ For 
example, the latter admit shock-like solutions as well as envelope light solitons 



Recently, motivated by observational and experimental studies |4|, |5|, |6|, B, Isj] , the governing 
equations for light interacting with electron-acoustic waves (EAWs) were derived j9|, using 
a two-species electron fluid model. Since ions are assumed immobile in Shukla et al. j^l, 
the phenomena of light-EAW interactions takes place on a timescale shorter than the ion 
plasma period. 

In this Brief Communication, we will develop a formalism suited for treating effects of 
partial coherence of the light source. The modulational instability properties of the resulting 
system of equations are analyzed, and the coherent and incoherent cases are contrasted. It 
is shown that partial coherence lead to damping of the modulational instability. However, 
in the limit of short wavelengths, the effect of a finite spectral width is suppressed, and 
such short wavelength modes are therefore more likely to yield, e.g. pulse filament at ion. 
The results can be of importance in space plasmas and laser-plasma systems involving short 
laser pulses. 

In Ref . |9( , the equations governing the nonlinear interaction between coherent light waves 
and electron-acoustic waves were obtained. We consider the propagation of electromagnetic 
waves in a two population electron plasma, where one of the electron species is hot, while 
the other electron species is cold. The electromagnetic wave is given by the vector potential 
A(t, z) = ip(t, z) exp(ikoz — iu t). The evolution of the slowly varying light wave envelope 
ip(t, z) is given by the equation 

(| - ^) * + c2 f^ - «&(JV fc + W = 0. (!) 

where v g = k c 2 /u is the group velocity of the light wave, N h = n h /n h0 , N c = n c /n c0 , 
<5 = nco/riho, n c (nh) is the cold (hot) electron number density perturbation, n c0 (riho) is 
the cold (hot) electron background number density, and cu p h = (Anrihoe 2 /me) 1 ' 2 is the hot 
electron plasma frequency. Moreover, the hot electron density distribution is given by the 
Boltzmann distribution, i.e. 

N h = exp(<^ - *) - 1, (2) 
where = ecp/Th is the normalized electron-acoustic wave potential and \& = e 2 \ip\ 2 /2m e c 2 Th 



m. 



is the normalized ponderomotive force potential due to the light source. The cold electrons 
are determined by 

8N 8 

^ + ^[(1 + ^K] = 0, (3a) 

and 

a +w "*h = ^- & ' (3b) 

where Vtk = {Th/m e ) l l 2 is the thermal speed of the hot electrons, and v c denotes the cold 
electron fluid velocity. Finally, the electron-acoustic potential is determined by Poisson's 
equation 

Akf^ = N h + N C S, (4) 

where Xdh— (Th/^rihQe 2 ) 1 / 2 is the hot electron Debye length. 

From Eqs. ©-© we obtain, using \tp — <C 1, the equation |t| 

where C e = VxhS 1 ^ 2 - Thus, Eqs. (JTJ), (jU) and (jSJ) form the desired system of equations for 
coherent interaction between light and EAWs. In Ref. it has been shown that (JTj), (J3J, 
and (jnj) admits localization and collapse of light pulses. 

In order to account for the effects of partial coherence in the light source, we employ the 
Wigner formalism. Starting from the two-point correlation function of the field of interest 
(in our case the light field), one performs a Fourier transform of this correlation function and 

n 

obtains a generalized distribution function [10(. Thus, we let the generalized distribution 
function for the light quanta be [ll]] 

p(t,z,k) = ±- J d(e ik Hr(t,z + (/2)^(t,z-(/2)), (6) 

where the angular brackets denotes the ensemble average. Then the light intensity is given 
by 

e 2 U\ 2 e 2 r 

Applying the time derivative on the definition (jBJ), and using Eq. ((TJ, we obtain the kinetic 
equation 

dp f c 2 k\dp <4A (d 2 p\ . (ldd. 
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for the light pseudo-particles. Here the arrows denotes the direction of operation, and the 
sin-operator is defined in terms of its Taylor expansion. By assuming a light distribution 
function with a finite spectral width, i.e. a finite spread in the light power spectrum, the 
effect of partial coherence can be incorporated into the light propagation through the two- 
electron plasma. Thus, Eqs. ©, (JZJ), and (jSJ) describes the interaction between partially 
coherent light and EAWs. 

Next we analyze the stability properties of the system ©, ((7j), and (jHJ) by perturbing 
these equations around a stationary plasma state (see, e.g. for a similar treatment of 
quantum plasmas). We let p(t, z, k) = po(k) + pi(k) exp(iKz — iflt), where \pi\ <C \p \, and 
ip = (fiexp(iKz — iQt). Moreover, we have \1/ = ^o + ^i exp(iKz — iQz), where j^xl ^ l^ol- 
Linearizing the system (J7J), and (jHJ) with respect to the first order quantities, we obtain 
the nonlinear dispersion relation 



2m e c 2 Th 2lj J ft — (v g + c 2 k/uJo)K 

where A = (ft 2 - C 2 e K 2 )/(tt 2 - C 2 K 2 + \ 2 Dh K 2 n 2 ). 

In the case of coherent light we have po = \ipo\ 2 5(k — k ), and the nonlinear dispersion 
relation Q yields 

t _ e 2 ujIA l^p (iQ) 



2m e c 2 T h 2c 2 [K -uj (tt-v g K)/c 2 K] 2 - K 2 /A 

i.e. 



(Q 2 -C 2 K 2 +X 2 Dh K 2 Q 2 ' 



ft- K + — )A 



2 4 7>-4 I ^2, ,2 \ 2 



4u 2 u 2u 2 



(11) 

Suppose now that the pulse phase &(z) experiences a random variation, which satisfies 
12| (exp[— i$(z + C/2)] exp[i$(2; — C/2)]) = exp(— A|£|), where 2A is the full wavenumber 



Po(*0 = ^ 71 ., XA2 , (12) 



width at half maximum of the power spectrum. Then the corresponding distribution function 
is given by the Lorentzian [l2[ 

W A 

7T (k - K ) 2 + A 2 

where k is a wavenumber shifting the location of the maxima of the distribution function 
Pok Inserting the expression (fT2j) into the nonlinear dispersion relation (JHl) we obtain 

1- ^ ^ A K 2 A- — — , (13) 



2m e c 2 T h 2c 2 [k - zA - cu (ft - v g K)/c 2 K} 2 - AT 2 /4 
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i.e. 



(Q 2 - C 2 e K 2 + X 2 Dh K 2 Q 2 



H - ( V g + — ^ " ] if 

CU 



n 2 



c 4 K 4 
4u 2 



■^ ".V ^n 8 - c 2 ^ 2 ), (14) 

2a; n 



Comparing the dispersion relations (Jllj) and (|14|) . it can be seen that the spectral broadening 

will reduce the growth rate. 

From now on we put the wavenumber shift kq to zero, thus centering the Lorentzian 

distribution (|T2j) around 0, as well as transforming to a comoving system, such that v g — > 0. 

We may use normalized and dimensionless variables defined by 

„ cK Q C e uj X Dh uj 2 h \ 2 Dh ^o C A 
K -> — , il -»• — , C e -> — , Ada -> , -> ^ , and A -> — . 

LUq UJq C C ZC z LUq 

(15) 

Using these dimensionless variables, we have plotted the typical behavior of the growth rate 
T = 2(ReO — Q), as given by the dispersion relations and (jl4j) . in Fig. ^ The uppermost 
curve corresponds to A = 0, i.e. the dispersion relation as given by (|10p. The growth rate 
asymptotically tends to 1 for the parameter values chosen in Fig. El and characteristically 
approaches an constant value for large wavenumbers if other parameter values are chosen. 
Moreover, the three lower curves has successively higher values of the spectral width A. The 
damping character of a nonzero spectral width can clearly be seen. Furthermore, it is clear 
that for long wavelengths, i.e. small values of K, the damping may completely suppress the 
growth rate, although small values may even enhance the growth rate (cf. the solid and 
dashed curved for small K). However, for short wavelength modes, i.e. higher values of K, 
even the cases with nonzero A asymptotically tends to the same value of the growth rate 
as in the monochromatic case. Thus, short wavelength modes seem almost unaffected by 
the spectral broadening, and such short wavelength perturbations would therefore dominate 
in the contribution to the filamentation, as well as in soliton and shock wave formation [l| 
of the incoherent light interacting with EAWs. In principle, these observations could be of 
importance in situations where the two-species electron model plays an important role, such 
as in space plasmas and laser-plasma systems , 0, 0, Q, 0] . 

To summarize, we have derived the governing equations for the interaction between light 
(coherent and incoherent) and EAWs. The limit of mono-chromatic light was analyzed, and 
compared to the case where spectral broadening was taken into account. It was shown that 



the effect of a finite width in the power spectrum of the light in general was to suppress 
the modulational instability. However, for short wavelength modes the spectral broadening 
was shown only to have a small influence on the modulational instability growth rate, thus 
making these modes more likely to dominate in the filamentation of pulses and the formation 
of shocks and solitons in space and laser plasmas. 
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FIG. 1: The growth rate V plotted as a function of the wavenumber K with the wavenumber shift 
Ko = 0. Here we have used the normalized variables (jl5|) . Moreover, we have used the values 
normalized values C e = 0.5, = 0.5, = 0.5. The solid black curve has A = 0, the dashed 
curve A = 0.1, the dotted curve A = 0.5, and the solid gray curve uses A = 1. The damping 
effect of a finite spectral width can clearly be seen, although the short wavelength modes are less 
affected. The asymptotic value of the growth rate in this case is 1. 
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